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ABSTRACT 

We measure the mean halo-mass of z ~ 0.5 Mgll absorbers using the cross-correlation 
(over co-moving scales 0.05-13/i^^ Mpc) between 1806 Mg II quasar absorption systems and 
~ 250,000 Luminous Red Galaxies (LRGs), both selected from the Sloan Digital Sky Survey 
Data Release 3. The Mg II systems have A2796 rest-frame equivalent widths ^^-^2796 > g.S A. 
From the ratio of the Mg II-LRG cross-correlation to the LRG-LRG auto-correlation, we 
find that the bias ratio between Mgll absorbers and LRGs is ^Mgll/^LRG = 0.65 ± 0.08, 
which implies that the absorber host-galaxies have a mean halo-mass ^ 20-40 times smaller 
than that of the LRGs; the Mgll absorbers have haloes of mean mass (logAfh(MQ)) = 
11.94 ± 0.31(stat)tQ25 (sys). We demonstrate that this statistical technique, which does not 
require any spectroscopic follow-up, does not suffer from contaminants such as stars or 
foreground and background galaxies. Finally, we find that the absorber halo-mass is anti- 
correlated with the equivalent width. If Mg II absorbers were virialized in galaxy haloes a 
positive Afh-W;:^^^^^ correlation would have been observed since ]/t/-'^2796 jg ^ direct measure 
of the velocity spread of the Mg II sub-components. Thus, our results demonstrate that the in- 
dividual clouds of a Mg II system are not virialized in the gaseous haloes of the host-galaxies. 
We review past results in the literature on the statistics of Mg II absorbers and find that they too 
require an My^-W^^'^^^ anti-correlation. When combined with measurements of the equiva- 
lent width distribution ((PN/dz/dWr), the Afh-VF/'^^^^ anti-correlation naturally explains 
why absorbers with VF/*^^^® ^ 2 A are not seen at large impact parameters. We interpret the 
Afh-W^:^^^^^ anti-correlation within the starburst scenario where strong Mg II absorbers are 
produced by supernovae-driven winds. 

Key words: cosmology: observations — galaxies: evolution — galaxies: haloes — quasars: 
absorption lines 



1 INTRODUCTION 

The connection between quasar (QSO) absorption line (QAL) sys- 
tems and galaxies is crucial to our understanding of galaxy evo- 
lution since QALs provide detailed information about the phys- 
ical conditions of galaxy haloes out to large impact parameters 
(p > 100 kpc) with no direct dependence on the host galaxy lu- 
minosity. Mg II absorbers are ideal for this purpose as the Mg IlAA 
2796,2803 doublet can be detected from z ~ 0.3 to z ~ 2.2 at 
optical wavelengths. Since the ionization potential of Mgl is less 
than 13.6 eV but the ionization potential of Mgll is greater than 
13.6 eV, Mgll absorbers trace cold gas. In fact Mgll absorbers 
with equivalent widths Wr^^^^^ > 0.03 A have been shown to 
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be associated with H I abs orbers covering five decades in H I col- 
umn density (A'j^j) (e.g. Petitiean & Bergeron 1990), including 
sub-Lyman limit s ystems I Churchill e t al. 1£9 9, 2000a), Lyman 
limit systems (e.g. iBergCTo n & Stasin sk31l98^[Steidel & S ar^enj 
'1992^ and damped Lya systems (DLAs; e.g.,'Le Bmi^^l^^9'^^ 
Rao & Turnshek 2000; Boisse et al. 1998; Churchill et al. 20001 
iRao et al J 12006) . which means that a large range of galactic en- 
vironments are likely to be sampled. For example, strong ab- 



with galaxies I'Lanzetta & Bowen"l99Cf; 'Bergeron & Boissa^^l 


Bergeron et al. 1992; Steidel & Sargent 


1992; Drinkwater et alj 


1993; Steidel etal. 1994). These groups 


lave shown that galaxies 



responsible for the Mgll host-galaxies have luminosities consis- 
tent with normal field 0.7 L% galaxies, and ISteidel et alJ showed 
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that, from their average colour, they are on a verage late-type (Sb) 
galaxies, a result reproduced bv lZibetti et alJ iioOS'). 

Medium-resolution spectroscopy of the Mgll absorbers 
quickly revealed that such absorbers are compos ed of sev- 
eral s ub-components ise rgeron & Stasin ska 198^ iTvtler et'all 
Il987t IPetitiean & BergerorJ Il990l) and that the number of 
sub-components strongly correlates with equivale nt vyidth 
fPetitiean & Bereeron 1990). Building on this work, Churchill 
U997..) showed that the mean Doppler width of individual compo- 
nents is 5kms^^, with an rms of comparable magnitude, using 
high resolution spectra (FWHM ~ 6kms~^). This co rresp onds 
to a thermal temperature of ~30,000 K. Churchill ( 19*3 also 
directly constra i ned th e turbulent component to be < 2kms~^. 
IChurchill et all <2003h . using high-resolution spectra, confirmed 
that equivalent width does str ongly correlate with the nu mber of 
sub-components as shown by IPetitiean & BergeroiJ fl990). This 
correlation arises because a large equivalent width can only be 
produced by more components spread over a large velocity range 
since extremely few components with large Doppler widths are 
seen. The equivalent width, W^^^^^^, should therefore be correlated 
with the velocity range, Au, covered by the sub-components and, 
indeed, this is observed to be the case (e.g. Ellison 2006, see 
their figure 3). The velocity range for strong Mg II systems are 
very large, from 50 to 400 kms~^. If the individual clouds were 
virialized within the haloes of the host galaxies. Aw would repre- 
sent the velocity dispersion of the many clouds in the host-galaxy 
halo. In this case. At; would be directly related to the mass of 
the host-galaxy. Since a larger W^'^'^^^ is achieved by having a 
larger number of clouds spread over a larger Av, W^.-^^Tse gj^^^^j^ 
also be positively correlated with the mass of the host galaxies if 
the individual clouds are virialized within the haloes of the host 
galaxies. 

Among the few correlations observed between the host- 
galaxy properties and the absorber properties. iLanzetta & BovJei] 
1 1990) noted a significant anti-correlation between Mgll equiva- 
le nt width and the impact pa ramete r distribution from the sam ple 
of lBergeron & BoissdiT99ll) . While lBergeron & Boissljl99ll) ar- 
gued that the ant i-correlation was not signi fican t, it was clearly seen 
in the sample of lSteidel & Sargend <I992h and lSteidell jl995l) : ab- 
sorbers with W,-^2796 > 2A are observed at small impact parame- 
ters but do not exist at large impact parameters (p > 50kpc). This 
is puzzling as large W^'^'^^^ absorbers at large impact parameters 
ought to be easily detectable. 

Only a slight correlation is observed between absorption 
cross- section and galaxy luminosity. ,Steid el l,1.993J) and Steidel 
i 1 9951) showed that the cross-section of Mgll absorbers with 

^A2796 ^ j^gg ^^(jjyg ^^^^^ ^ 4Q /j-l]^pc (physical, 

70 kpc co-moving) from the observed absorber-host-galaxy 
impact parameter distribution. Furthermore, they concluded that 
the cross-section radius only slightly increases with luminosity, 

-Rphys CC {Lk/L*j() ■ . 

The only other known correlation between galaxy proper- 
ties and absorber properties is that recently discovered between 
^/A2796 and the galaxy morp hological asymmetry i Kacprzak et alj 
l2005ll200d) . lKacprzak et alJ searched for correlations between gas 
properties, host galaxy impact parameters, inclinations (& position 
angles) and morphological parameters. Among those parameters, 
they only found a 3.2-a correlation between W^^"^^^ and the asym- 
metry in the host galaxy morphology as measured from the residu- 
als of 2D light profile fits. 

From studies of individual absorbers, some recent work that 
compared the absorption kinematics with the host galaxy rotation 



curve favour the idea that strong Mg II cloud s arise in galactic out- 
flows jBond et alJl200ll : lEUison et alj|2003h . On the other hand, 
'Stei del et alj i2002') showed that, in 4 out of 5 Mgll absorbers 
(selected to be edge-on and aligned towards to QSO), the absorp- 
tion kinematics are consistent with rotation being dominant for the 
absorbing-gas kinematics. However, a simple extension of the rota- 
tion curve fails. 

Many authors h ave tried to put QALs in the context of the- 
oretical models since lBahcall & Spitzer 1 1969) who proposed that 
the metal absorption lines are produced by gaseous haloes of inter- 
vening galaxies with a large cross-section (up to 100 kpc). Later, 
[York et al. 1 1986) argued that QA Ls arise mainly in the haloes 
of gas rich Magellanic -type dwar fs. Mo & Miral da-Escudl il996ft 
followed the ideas of lBahcall & Spitzer ll969i) and produced a 
detailed model in which Mgll absorbers are signatures of in- 
falling (photo-ionized) cold gas embedded in a T ~ lO" K gas 
halo. In this model, the in-falling cold gas should be virialized. 
[Mailer &Bullock|j2004l) reached similar conclusions. 

Despite these numerous past results, a fundamental question 
remains: What is the physical nature of strong Mgll absorbers? 
Here we constrain one important physical property of Mg II ab- 
sorbers - namely the halo-mass of the host galaxies - statistically 
by studying the clustering of galaxies around the absorbers. Clus- 
tering studies of metal-line systems, such as the absorber-absorber 
auto-c orrelation in velocity space, have been used for some time 
(e.g. Sargent et al. 1988';'Steidel & Sargent ' 1992; Churchill ' 1 997|; 
.Charlto n & Churchill 1998). This technique measures the line-of- 
sight clustering and therefore suffers from strong peculiar veloc- 
ity effects. Numerical models are then required to infer the physi- 
cal properties (halo mas s, halo sizes, etc.) of the host galaxies (see 
IScannapieco et a l. 2006. for a recent example). 

In order to avoid these limitations, we choose to cross- 
correlate Mg II absorbers with 'field' gala xies tracing the two- 
dimensional large scale structure as in Bouche & Lowenthp 
^004), Bouche. Murphy & Peroux (2004), Bouche et al. (200^ 
ICookeetal . 12006) and lRvan- We ber (2006). Specifically, we mea- 
sure the amplitude ratio of the absorber-galaxy cross-correlation 
to the galaxy-galaxy auto-correlation. In hierarchical galaxy for- 
mation scenarios this is a direct measure of the ratio between the 
bias of the absorber host-galaxies and of the field galaxies used. 
From this bias ratio, the mass of the absorber host-galaxies can be 
inferred. 

In this paper, we extend our DRl results iBouche et all2004l) 
by cross-correlating 1806 Mgll absorbers with ~ 250,000 Lumi- 
nous Red Galaxies (LRGs), both selected from the Sloan Digital 
Sky Survey (SDSS) Data Release 3 (DR3; Abazaiian et al. 200f|), 
over ~ 800 square degrees. Thanks to the large SDSS survey, this is 
a leap for ward from other clus tering studies such as Williger et al] 
(2002) and lHaines et a l ( 2004) which cover a few square degrees. 

This paper is organized as follows. In Section|2| we summarise 
the selection of our Mgll absorbers and LRGs. In Section |3| we 
describe our method to measure the halo-mass. The results are pre- 
sented in Section |4] We test these results against numerous past 
results on Mgll absorbers in Section|5|and discuss a physical in- 
terpretation of our results in Section |6| Our main results are sum- 
marised in SectionQ For those familiar with absorber-galaxy clus- 
tering analyses, a quick read of this paper comprises Fig. [S] and 
Fig.jS] followed by SectionjS] A critical reader should focus on the 
several consistency tests we performed (Figs. |4| & |6| and on the 
discussion in SectionjS] 

We adopt D,m = 0.3, Ha = 0.7 and Hq = 
WOh kms~^ Mpc~^ throughout. Thus, at z — 0.6, 1" corresponds 



Strong zc::^ 0.5 Mg II absorbers 3 



to 7.44/i ^ kpc and (52 = 0.1 corresponds to 216/i ^Mpc, bothin 
co-moving coordinates. 



2 SAMPLE DEFINITIONS 
2.1 Mg II absorbers 

The algorithm to select Mgll absorbers from SDSS/DR3 dif- 
fers in several way s from the algorithm used for SDSS/DRl in 
iBouche et alJ i2004^. the most important of which is the method for 
estimating the QSO continuum. In Bouche et al. (2004) we used a 
series of overlapping polynomial fits to small Sections of the con- 
tinuum. While this provides a reliable continuum in most cases, it 
does not perform well near sharp QSO emission lines, particularly 
when absorption lines - possibly the target Mg II lines - are im- 
printed over the emission. To alleviate this problem in the current 
SDSS/DR3 analysis, we used principal component analysis (PCA) 
reconstructions of the QSO continua. Full details of the method are 
reported in Wild et al. ( 2006, see their appendix). 

The features of the PCA algorithm most important for the 
present paper follow. Eigenspectra were generated in 4 QSO emis- 
sion redshift bins to reduce the amount of data 'missing' due to 
the differing wavelength coverage of each spectrum: 0.005-0.458, 
0.381-0.923, 0.822-1.537, 1.410-2.179, 2.172-3.193. An iterative 
procedure was used to identify and remove quasars showing broad 
absorption lines (BALs) during the creation of the PCA eigen- 
spectra. This improves the continuum reconstruction of non-BAL 
quasar spectra, by removing from the input sample features which 
vary greatly in a small number of objects. Fig.Qshows an example 
SDSS spectrum with its PCA continuum. 

Having established a continuum, candidate Mgll absorbers 
are identified as follows. We searched for intervening Mgll ab- 
sorbers from Zabs = 0.37 to z^bs ~ 0.8 (see Section IZ2l . The 
low redshift cut arises from the fact that SDSS QSO spectra be- 
gin at ~3800A, and the high redshift cut was imposed since this is 
where LRGs drop below the SDSS magnitude limits. In addition, 
we require all Mg II systems to be above the Lya QSO emission 
line. QSOs above Zcm ~ 3.193 are therefore not considered. All 
pixels above 1250 A in the QSO rest-frame and > 3000 km s^^ 
blue-wards of the Mgll emission line are tested for Mg II A2796 
absorption. At each pixel, putative Mg II A2796 lines are charac- 
terized using a method similar to that detailed by Schneider et al. 
( 1993). For initial line detection, the rest-frame Mg II A2796 equiv- 
alent width and l-a detection limit are defined using the spectro- 
graph instrumental profile (IP) as a weighting function and using 
pixels in a 7 A window in the putative absorber's rest-frame. The 
IP is assumed to be a Gaussian of width FWHM — 160 kms~^. 
If this estimate of the Mgll A2796 equivalent width is ^ 0.3 A 
and is significant at ^ 8 cr then the putative absorption redshift is 
estimated from a parabolic interpolation of the equivalent widths 
of the current pixel and its 2 neighbouring pixels. The equivalent 
width of the Mg II A2803 line is estimated in a similar way based 
on the A2796 redshift. If the A2803 equivalent width is significant 
at ^ 3 (T then the system is flagged as a candidate Mg II absorber. 
Spurious candidates are removed by visually inspecting each Mg II 
candidate. The most common mis-identification is broad C IV ab- 
sorption near the C IV emission line. Adopting a conservative ap- 
proach, we rejected any candidates which did not show absorption 
in either Fe II A2600 or Mg I A2852 at ^ 1 cr significance. 

For each candidate we derive a refined estimate of the Mg II 
A2796 equivalent width using a Gaussian fit to the absorption as a 



Table 1. Catalogue of 1806 Mgll absorbers from the SDSS DR3 with 
0.37 ^ 2abs ^ 0.80. The J2000 name, QSO and absorption redshifts 
and the measured Wr for Mg IIAA2796 & 2803, Mg I A2852 and Fe II 
A2600 are given. Here we show only a small sample from the full 
table which is available in the electronic edition of this paper and from 
'http://www.ast. cam. ac.uk/^mim/pub.html Full name designations, and 
statistical uncertainties in Wr are given in the electronic version, together 
with Wr measurements for the following transitions: Zn II AA2026/2062, 
Cm A2056/A2062/A2066, Fell A2344/A2374/A2382/A2586, Mnll 
A2576/A2594/A2606, Till A3242/A3384, Call AA3934/3969 and Nal 
AA5891/5897. 



Rest equivalent width 



SDSSJ 




^abs 


2796 


2803 


2852 


2600 


160530+393116 


1.083 


0.4969 


1.04 


0.79 


0.04 


0.19 


160726+471251 


1.816 


0.4974 


1.30 


1.41 


0.03 


0.59 


171726+654542 


1.215 


0.4974 


3.29 


3.39 


0.22 


1.39 


112719+654143 


1.250 


0.4977 


0.53 


0.32 


0.17 


0.24 



weighting function (c.f. the IP-weighting above). This new, some- 
what more optimal, estimate is referred to throughout this paper as 
the measured Mg II A2796 equivalent width, W^^'^'^^^, for the sys- 
tem. Equivalent widths, with a similar Gaussian-fit weighting, are 
also derived for a variety of other commonly observed transitions 
(see Table0. 

With the above algorithm we detected and visually confirmed 
1806 Mgll absorbers in SDSS/DR3. Figure IH (left panel) shows 
the distribution of the equivalent width of our Mg II absorbers. Ta- 
bleQis an excerpt from the catalogue of absorbers which is avail- 
able in its entirety in the electronic edition of this paper and from 
an on-line catalogue at http://www.ast.cam.ac.uk/~mim/pub.htrn^ 
FigureQshows an example Mg II absorption system. Similar plots 
are available for all absorption systems in the on-line catalogue. 

2.2 Luminous red galaxies 

lEisenstein et alj I 200 ih (see also IScranton et al."2003') presented 
colour criteria specifically designed to select luminous massive 
early types both locally and at z > 0.4, i.e. much beyond the vol- 
ume of the SDSS Main sample. Because of their association with 
luminous {Air ~ —21.5), massive haloes and spectral uniformity, 
LRGs are excellent probes of the large sc ale structure as proven b y 
the detection of the baryon oscillations bv lEisenstein et alii 2005^ . 

For each Mg II absorber, galax ies meeting the following crite- 
ria (foUowing'Scrant on et all2003h were extracted from the SDSS 



DR3 galaxy catalogue: 

*pctro < 21, (1) 

0.7(g* -r*) + 1.2[(r* -r)-0.18] > 1.6, (2) 

{g'-r') > 1, (3) 

(dx EE)(r* -i*) - (3* -r*)/8 > 0.4, (4) 

' model 

> 0.24 , (5) 

IZphot— Zabsl < 0.05. (6) 



We also required errors on the model magnitudes to be less than 
0.2 mag in r* and i*, and we excluded objects flagged by SDSS 
as BRIGHT, SATURATED, MAYBE.CR or EDGE. The model 
magnitudes, corrected for Galactic extinction, were used to com- 
pute the colours. Equations Q-JIJ ars the LRG selection crite- 
ria of Scranton et al. (2003). Criterion ^4} is equivalent to impos- 
ing 2phot 0.3. Criterion Js} separates stars from galaxies. Crite- 
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Figure 1. SDSS spectrum of QSO SDSSJ02521 1.9 1-080244.2 containing a strong Mgll absorber at Zabs = 0.7688. The upper panel shows the entire 
spectrum (black histogram) with our PCA reconstruction of the QSO continuum (solid line). The error array is shown by the continuous line just above the 
zero flux level. The vertical dashed lines mark several transitions at the absorption redshift, including the Fe II A2600 and the Mg I A2852 lines that we 
used to confirm the Mgll doublet. These transitions are detailed in the lower panels. Detailed plots for all of our 1806 Mgll absorbers can be viewed at 
http://www.ast.cam.ac.uk/~mim/pub.html. 
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Figure 2. Left: Rest-frame equivalent width (H//'^^^^) distribution of our 1806 Mgll systems. Right: Redshift distribution of the LRGs (bottom) and Mgll 
systems (top). The inset shows the distribution of QSO emission redshifts. 



rion j6| is the selection of galaxies within a redshift slice of width 
Az =0.1 around Zabs using the photometric redshifts, Zphot, of 
ICsabai et aljj2003h who showed that these are accurate to = 0.1. 

The choice of the slice width =0.1 corresponds to ~ 
200/i"^ Mpc (CO -moving) and is arbitrary. We will show that the 
amplitude ratio between the absorber-galaxy cross-correlation and 
the galaxy-galaxy auto-correlation does not depend on the choice 
of Az in Section l4!4l since we use the same redshift width for both 
correlation functions. 

Finally, we remove the '--^ 10 per cent of galaxies with prob- 
lematic photometric redshifts by requiring that galaxies have Zphot 
uncertainties CTz < 0.5. A total of 242,620 galaxies met all these 
criteria in our 1806 fields (~ 800 square degrees). 

Figure fright) shows the redshift distribution of these LRGs 
for the 1806 fields. We used the spectroscopic redshift when avail- 



able, which includes less than 1% of the sample. This small fraction 
is likely to increase in future with the joint 2dF/SDSS program to 
obtain spectra of LRGs (e.g. Padmanabhan et al. 2005) . 



3 METHOD 

We first describe the basics of the galaxy clustering analysis in Sec- 
tion The correlation estimator best used for this work is dis- 
cussed in Section lT2l 

3.1 Galaxy clustering around QSO absorbers 

A widely used statistic to measure the clustering of galaxies is the 
correlation function, ^(r). The absorber-galaxy cross-correlation. 
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^ag, is defined from the conditional probability of finding a galaxy 
in a volume dV at a distance r — |r2 — ri|, given that there is a 
Mg II absorber at ri : 



P(LRG|Mgii) = n„ [1 + eag(r)] dV, 



(7) 



where n„ is the unconditional background galaxy density. 

Because the observed amplitude of the auto- and cross- 
correlation functions are related to the dark matter correlation func- 
tion, 5dm, through the mea n bias, b(M), which is a function o f the 
dark matter halo-mass fe.g. lMo et alll993l:lMo & Whitel20ol and 
references therein). 



Cgg(r) = 6(MgKDM(r), 
5ag(r) = 6(Ma)6(Mg)eDM(r). 



(8) 
(9) 



the relative amplitude of the cross-correlation will give a measure- 
ment of the bias ratio and thus of the relative halo-masses: 



6(Ma 



(10) 



In other words, if the amplitude of ^ag is greater (smaller) than ^gg, 
the haloes of the absorbers are more (less) massive than those of 
the LRGs since the bias increases with halo mass in all hierarchical 
models. 

In the remainder of this work, we will use only projected 
correlation functions, w{rg), where re — Da(1 + z)9 for Da 
the angular diameter distance in co-moving Mpc. This is neces- 
sary since our sample is made up of absorbers with spectroscopic 
redshifts and of galaxies with photometric redshifts. The projected 
cross-correlation between Mgll absorbers and LRGs, uiag(re), is 
related to Cag(r) via Wag(re) = J A^(z)Cag(\A^ + ^) dl where 
N{z) = dN /dl is the line-of-sight distribution of LRGs and I is 
the co-moving distance along the line of sight. 

For galaxies distributed along the line of sight as a top-hat 
function of width (normalized such that J dN/dl dl = 1), the 
amplitudes of both tUgg (rg) and Wag (rg) are inversely proportional 
to Az (see also JEisenstein 2003; Bouche et al. 2005, appendix A): 

"'gg(''e) = ^gg '■e"^ - ^e"^ '^o.gg ^7 X ^ > (H) 



5 (re) = yla 



1 — / 1 — / 7 



(12) 



where 



r(*)r(^) 



— 'Fr^\ ^f"l '"o.gg & '"o,ag are, respectively, the 

galaxy-galaxy & absorber-galaxy correlation lengths. These two 
equations show that both uiag and Wgg depend in exactly the same 
way o n the width of the redshift distribution Az . Bouche et al. 
y005, appendix A) showed that this is always true when one cor- 
relates one population with known redshifts (the Mg II absorbers) 
with another population whose redshift is more uncertain. In Sec- 
tion l4.4l we will show empirically that the ratio Wa,g/wgg is indeed 
independent of A^. Thus, we stress that Wn^/w^g is independent 
of the width of the LRG redshift distribution as long as one uses 
the same galaxies for lUag and for w^g . 

From equations <1 1> & <12> . the ratio of the amplitudes of the 
two projected correlation functions. 



Wa 
We, 



(13) 



is simply (ro.ag/ro.gg)'', and is also equal to the bias-ratio 
&(M]y[g jj)/6(Mlrg) [equation jlOU from which we infer the 
mean Mg II halo-mass using the bias prescription of IMo & Whitd 
j2002h. 



It is important to realise that measuring the halo-mass from 
the ratio of projected cor relation functions has the following ad- 
vantages, as advocated in iBouche et alj jjOoT) and Bouc he et alJ 
(2005): (i) one constrains the mass of the Mgll host-galaxies in a 
statistical manner without directly identifying them, (ii) it is free 
of possible systematic errors due to foreground or background con- 
taminants (i.e. stars or galaxies) and (iii) it does not require knowl- 
edge of the true width of the redshift distribution of the galaxy 
population. The first point is a natural consequence of correlation 
statistics. The last two points are consequences of the fact that the 
same galaxies are used to calculate 'Wgg{rg) and Wa,g{rg) and that 
both Wag and w^g have the same dependence on the width of the 
galaxy redshift distribution [equations <1 H & <12H . This last point 
follo ws from the fact tha t the absorber redshift is known precisely 
(see iBouche et all200^ appendix A) as mentioned above. We will 
demonstrate points (ii) & (iii) empirically in Section l4!4l 

Usin g smoothed particle hydrodynamical (SPH) simulations 
IBouche et alJ ^20051) showed that the measured amplitude ratio 
liiag/uigg returns the bias-ratio ba/bg that one expects given the 
known m ass of the simulated galaxies. Specifically, the bias for- 
malism o f lMo&Whit3 i2002') predicts an amplitude ratio of 0.771 
in the case of the DLA-galaxy cross-correlation - a situation 
very similar to the Mg Il-galaxy correlation considered here. Di- 
rect measurement of the correlation functions from the simulations 
yields boLA/bgai = 0.73 ± 0.08, in excellent agreement with 
the prediction. The simulations also demonstrated that the cross- 
correlation technique does not depend on the galaxy population 
used to trace the large scale structure and that foreground and back- 
ground contaminants do not effect the measured bias-ratio. 



3.2 Which correlation estimator to use? 

Which estimator is best used to compute the Mgll-LRG cross- 
correlation? For a given field (i.e. a single absorber), a seemingly 
natural choice would be the estimator 

1 + t^ag [rg ) = 



R 



(14) 



where G is the observed number of galaxies between rg — dr/2 and 
rg + dr/2 around the absorber, and R is the number of absorber- 
random galaxy pairs normalized to the total number of galaxies in 
the field, A'^g, i.e. multiplied by Ng/Nr. Ideally, R would simply 
be the area of the annulus times the surface density of galaxies. 
In practice, some fields are on the edges of the SDSS coverage, 
and the SDSS coverage itself is not completely uniform on the sky; 
there are small holes, areas missed by SDSS or not yet released. 
One overcomes these problems by (1) generating random galaxies 
excluding the gaps and edges, and (2) by generating ~ 200 times 
more random galaxies than real galaxies in order to reduce the shot 
noise in R to an insignificant proportion. 

In our case, we have 1806 fields, at relatively precise absorp- 
tion redshifts, ^abs- The selection function for the fields is therefore 
a sum of (5-functions. Should one then use (■^), where () represents 
the averaging operator, in the right-hand side of equation J14t ? This 
is clearly not optimum since it would treat eac h field equally in 
performing the average. 'Adelbe rger et alj (|2003) (see their equa- 
tion B3) showed that a better choice for the estimator of u;ag(rfl) is 
{G)/{R),i.e.: 

AC 

l + wag(re) = — , (15) 

where AG — ^ , Gi is the observed number of absorber-galaxy 
pairs between rg — cir/2 and rg + dr/2, summed over all the fields 
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i. AR is the normalized number of absorber-random galaxy pairs. 
The normalization is applied to each field independently: AR = 
y^ ^. R* Nl/Nr, where R' is the number of random pairs in field i, 
and Ng (7V^) is the total number of galaxies (random galaxies) for 
that field. Naturally, we took into account the areas missing from 
the survey within our search radius. Note each absorber's redshift is 
used to compute re from a pair separated by an angular separation 

e. 

We stress that the lLandv & Szala^ il993h estimator is not ap- 
plicable here, since it is symmetric under the exchange of the two 
populations by one another. Here, the absorbers have spectroscopic 
redshifts, whereas the galaxies have photometric redshifts breaking 
the symmetry. Note that this asymmetry has an important conse- 
quence: it implies that both xuag and Wgg have the same dependence 
on the gala xy redshift distribut ion, as indicated by equations Jilt 
& (12} (see lBouche et all2005L appendix A). 

The cross-correlation computed from equation <15t is biased 
low due to the integral constraint (e.g. Peebles 1980), explained be- 
low. The reason for this is simply that the true correlation is defined 
as an overdensity with respect to the 'unconditional' galaxy density 
n„. This is not a measurable quantity as one uses the 'observed' 
galaxy density (or surface density Eg for projected correlations) 
for the field over which one measures w{re). As a conseqence, the 
sum of all the pairs must be equal to the total number of galaxies Ng 
observed in the survey of area E. This implies that the integral of 
uiag over the area E vanishes: Wa,g dE — 0. To fulfill this con- 
dition, the cross-correlation will have to be < on the largest 
scales, i.e. biased low. This bias is refered to as the 'integral con- 
straint', C. The true correlation function is then Wag = ifag + C, 
where C is the integral constraint. 

There are two ways to estimate C. Firstly, it can be estimated 
iteratively. For a given amplitude of the correlation, one can calcu- 
late the expected bias C which is then used to correct the estimate 
of the amplitude. The second method is to fit C when one has suf- 
ficient signal to noise. In the analysis below, we find that the first 
method gives 0.0217, while the second gives C = 0.018 ± 0.012 
(keeping the other parameters fixed). 



4 RESULTS 

In this Section, we first present the Mgll-LRG cross-correlation 
(Section 14. 1> . and the LRG-LRG auto-correlation (Section 14. 2t 
before showing the main result on the amplitude ratio between 
the Mg II-LRG cross-correlation and LRG-LRG auto-correlation 
(Section f4.4> . We turn the amplitude ratio into a halo-mass for the 
Mg II absorbers in Section l431 Finally, we show that the amplitude 
ratio (or the halo-mass) varies with equivalent width in Section l431 



4.1 Mg II-LRG cross-correlation 

Figure|3|(filled circles) shows the Mg II-LRG cross-correlation for 
the entire sample of 1806 Mgll absorbers, where we used equa- 
tion <15> for the estimator of Wag(r'e). There are 242,620 objects 
within rg — 12.8/i~^ Mpc which is the outer radius of the largest 
bin used. Table|2|shows the total number of pairs, AG, and the ex- 
pected number of pairs, AR, if Mg II absorbers and LRGs were not 
correlated. Figure |3] demonstrates that Wa,g{re) is a power-law at 
all scales. As a consequence, the d ip at ^ 100 kpc in the DRl 
sample discussed in Bouche et^Jiooi) appears to have been due 
to small number statistics. 
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Figure 3. The filled circles show the projected Mg II-LRG cross- 
correlation, Ws^g{rg), for 1806 absorbers and ^ 250,000 LRGs and the tri- 
angles show the LRG-LRG auto-con'elation, Wgg(rg). The error bars are 
computed using the Jackknife technique. The independent fits to uiag {rg ) 
(solid line) and Wgg{rg) (dashed line) are shown. We fixed the integral 
constraint C to 0.02 (see text). 

Table 2. The total number of absorber-LRG pairs, AG, and the number of 
absorber-random pairs, AR, for the cross-correlation shown in Fig.l3l 



rg AG AR uiag cr^ 

[h-l Mpc] 



0.05-0.1 


28 


11.61 


141 


0.42 


0.1-0.2 


82 


45.13 


0.82 


0.19 


0.2-0.4 


262 


181.92 


0.44 


0.10 


0.4-0.8 


900 


724.91 


0.24 


0.044 


0.8-1.6 


3183 


2904.27 


0.096 


0.015 


1.6-3.2 


12205 


11619.1 


0.050 


0.010 


3.2-6.4 


46867 


46043.9 


0.018 


0.007 


6.4-12.8 


179093 


181089 


-O.OIl 


0.002 



T he error bars for Wag are computed using the jackknife esti- 
mator <EfroiJll982h : we divide the sample into 10 parts and com- 
pute the covariance matrix from the A'jack = 10 realisations for 
each part: 

GOV,, = " ^ V [Wk{re^)-w{rg^)]-[wk{rg^)-w{re^)](l6) 

where Wk is the fcth measurement of the cross-correlation and w is 
the average of the A'jack measurements. 

We fitted Wagire) with a power law model, 'Wa_g{rg) = 
^ag J'e , by minimizing 

x' = -^[w-w]^COV-i[w-w], (17) 
Adof 

where A'dof is the number of degrees-of-freedom, w and w are 
the vector data and model respectively, and COV~^ is the inverse 
of the covariance matrix. Since GOV is singular we used singu- 
lar value decomposition techniques to avoid instabilities in its in- 
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Figure 4. The Mgll-LRG cross-correlation amplitude as a function of 
the artificial redshift offset, ^offset, iinposed on the absorber redshifts. 
As expected, the cross-correlation signal vanishes when l^offsctl is if" 
creased. This shows that the cross-correlation signal in Fig. |3] is solely 
due to the Mg II absorbers. The solid line represents a top-hat redshift 
distribution of width Az = 0.1 convolved with a Lorentzian of width 
FWHM = 0.176, corresponding to the typical photometric redshift un- 
certainty of cTz = 0.075. The solid curve is scaled to the ^offset = data 
point. The dotted line shows the sohd line shifted by Zshift = —0.025, rep- 
resenting the possible 5 per cent corre ction to the photometric redshifts sug- 
gested bvfPadmanabh anet alj |2003) (see text). The increase in the eiTor- 
bars from left to right is due to the shape of the LRG redshift distribution in 
Fig.llright). 

version (see discussion in iBemsteirl 1994) . Since the integral con- 
straint is C = 0.02 ±0.01 (see last Section), we add (5^^ (0.01)^ to 
the covariance matrix COV(w)ij to form COV(it)'")i,j, the co- 
variance matrix for {re) = w{rg) + C. 

Fitting the vector w"^, the best-fit amplitude at lh~^ Mpc and 
power-law slope of the cross-correlation are, respectively, 

y4ag = 0.143 ±0.012, 
/?ag = -0.893 ± 0.043. 



4.2 Luminous red galaxy auto-correlation 

The LRG-LRG auto-correlation «;gg is represented by the grey/red 
triangles in Fig.|3| Since our goal is to measure the ratio Wag/wgg, 
we are forced to use the same estimator for tUgg as lUag, namely 

. , GG ^ , „^ 

l + w^^{rg)^ (18) 

where GG is the observed number of galaxy-galaxy pairs between 
rg — dr/2 and rg + dr/2 and GR is the number of galaxy- 
random galaxy pairs, computed as before. Again, rg is in units of 
co-moving Mpc. The errors and the covariance matrix for Wgg are 
computed using iVjack ~ 10 jackknife realisations, as before. Using 
the covariance matrix to w^^, COY {w^)ij (as in Section 4.1), we 
fitted a power-law, Wgg{rg) = Agg r^*^*, to 'Wgg{rg). The best-fit 



amplitude at Ih ^ Mpc and the power-law slope are, respectively, 
Agg = 0.186 ±0.012, 

/3gg = -1.023 ±0.035. (19) 

The conversion of the amplitude, Agg, to the co-moving cor- 
relation length, r-Q.gg, depends on having precise knowledge of the 
true redshift distribution, N{z), of the LRG sample. The observed 
distribution is a convolution of the true distribution and the pho- 
tometric redshift errors. Deconvolving these to find the true dis- 
tribution is difficult (e.g. Padinanabhan et al. 2005) and is beyond 
the goal of this paper. Howev er, we note that ove r a redshift range 
similar to ours (z = 0.3-0.9). iBrown et ai]j2003h showed that red 
{Bw — R> 1.44)^ galaxies in the NOAO deep wide survey have 
a correlation length of ro = 6.3 ± 0.5h~^ Mpc, for galaxies in 
the luminosity range — 21.5 < Mr < —20.5. At the mean redshift 
of our sample, z = 0.5, such clustering is consiste nt with halo- 
masses of 1 x 10*^^ M p, using the bias prescription o f :Mo & White! 
( 2Q02j). lBrown et all showed that the correlation length rapidly in- 
creases to r'o,gg = 11.2h~^ Mpc at Mr — —22. Our sample has 
a mean luminosity of Mr* — —21.5, which is consistent with 
T'o.gg = 8h~^ Mpc, and a halo-mass slightly higher: Mlrg ~ 2- 
4 X 10^^ Mq . Since we will see that the systematic errors in our 
final results are of the same order as the statistical errors, we here- 
after assume^ that Mlrg ~ 3x 10^^ M0 and treat the uncertainty 
in its value as an additional source of systematic errors. This is dis- 
cussed further in Section l431 

4.3 Is the cross-correlation signal really due to the 
absorbers? 

In order to verify that the signal of the cross-correlation is solely 
due to the Mgll absorbers, we repeat our cross-correlation mea- 
surement adding an artificial offset Zofisot to the absorber redshift 
-Zabs, with ^offset ranging from —0.25 to ±0.25. Here, we fit the 
cross-correlation liiag(r'e) using the model 'w{rg) — A (r*^ — 
0.02/0.15) in order to account for the fact that as A tends to 
zero, so does the integral constraint C. The slope l3 is kept fixed 
at/3 = -0.89. 

Figure |4] shows the amplitude of the Mg II-LRG cross- 
correlation Aag as a function of the artificial redshift offset ^offset ■ 
Because the amplitude vanishes when |zoffsct | is increased, this fig- 
ur e shows that the Mg I I-LRG cross-correlation signal in Fig.|3|and 
in lBouche e t al.' ('2004') is solely due to the Mg II absorbers. 

Given that Csabai et al. |2003) showed that the LRG photo- 
metric redshifts have a typical uncertainty of = 0.05-0.1, we 
convolved our top-hat redshift selection function of width Az — 
0.1 with a Lorentzian with FWHM = 2.35 cr^, representing the 
typical photometric redshift uncertainty with an underlying popula- 
tion of 'outliers' . The solid line in Fig.|4]shows the result of the con- 
volution with az — 0.075. The data points in this figure might in- 
dicate that the photometric redshifts are slightly over-estimated. In 
fact, .Padmanabhan et al. 1 2005) showed that at redshifts z > 0.4, 

^ This cut is similar to ours, i.e. g — r > 1.0: According to the photometric 
transformations listed on the SDSS website, Bw — i? > 1.44 corresponds 
toB -V > 1.24.01 g- r > 1.14. 

^ After completion of this analysis and submission of this paper, 
[Mandelbauin et al. ( 2006) presented a tensing measurement of the shape of 
the density profile of galaxy groups and clusters traced by 2: < 0.3 LRGs. 
Their results imply that LRGs fainter than Mr = —22.6 (corresponding to 
our sample) reside in haloes of mass (2.9 ±0.4) X lO^^ M0. 
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Figure S. Same as in Fig.|3]in the two-halo regime rg > 1 Mpc. The 
dashed line represents «)gg and the relative amplitude of tOag to uigg (found 
using iDag = a X i«gg) is a = 0.808 ± 0.096. From this, the bias-ratio 
''Mgll/^LRG is ^ 0-65 (see text), which in turn implies that Mgll host- 
galaxies have haloes of mass A/h ~ 10^2 (see SectionP^ 



photometric redshifts are slightly over-estimated by 5 per cent, i.e. a 
redshift correction of 0.025 at the mean redshift of our sample. The 
dotted line shows the shifted curve. We emphasize that this will not 
affect our results since we measure Wag/^J^gg, which has no depen- 
dence on the redshift distribution of the galaxies. 

The error bars in Fig. |4| increase strongly with increasing 
-^offset • This is due to our redshift selection function shown in Fig.|2| 
as a positive offset is added, far fewer galaxies are selected and the 
signal-to-noise ratio decreases since the number of galaxies directly 
drives the number of absorber-galaxy pairs. 



4.4 The relative amplitude lOag/iOgg 

In order to measure the Mg II-LRG bias-ratio [from equations (S) 
& j9j], one needs to use (i) the scales where the bias b{M) dom- 
inates (i.e. where the correlation arises from 2 different haloes) 
and (ii) the same power-law slope, /3, in order to compare the 
amplitudes of Wag and Wgg. The first point requires that we use 
only scales re > lh~^ Mpc, as numerous papers (both on sim- 
ulations and on SDSS data) have shown that the correlations be- 
low these scales are dominated by the single- halo correlation be- 
tween the central galaxy and its satellites (e.g. lBerlind et all2003t 
[Zehavi et al. 2004). The auto-correlation of z = SDSS galax- 
ies shows a break at Ih^^ Mpc where the transition between the 
single- and the two-halo terms occur. Point (ii) is easily achieved 
by using Wgg (Section l4.2t as a template to constrain the relative 
amplitude of Wag : 



a X T^gg, 



(20) 



where a is the amplitude ratio. 

Figure|5|shows the auto- and cross-correlation on scales larger 
than re > 1 Mpc. We find that the best relative amplitude is 




0.1 



0.15 



0.2 



Figure 6. Top: Amplitude ratio a [equivalently, the bias-ratio 
b{M!i)/b{Mg)] as a function of the width Az of the redshift slice. 
The bias-ratio is independent of the redshift width selection. As a conse- 
quence, the ratio of the amplitude of the cross-correlation to that of the 
auto-correlation is free of systematics from contaminants (stars or inter- 
loping galaxies). Bottom: The amplitude Aag of the cross-correlation uiag 
decreases as the redshift slice selection increases, as expected. The vertical 
dashed line shows = 0.1, the selection adopted here [equation l6ll. 
Both panels were computed in the rg > 1 Mpc regime as in Fig.l3l 



: 0.808 ± 0.096 . 



(21) 



As we emphasized in lBouche et alJ <2004 and pointed out in 
Section l4m the same galaxies are used to calculate w^g and Wgg. 
Therefore, the relative amplitude, a, is free of systematics from 
contaminants (stars or interloping galaxies). This is demonstrated 
in Fig.|5| The top panel shows a as a function of the width of the 
redshift slice. As one increases up to 0.2 (~ 400 Mpc), i.e. as 
one increases the number of foreground and background galaxies, 
the amplitude ratio is independent of that choice. The bottom panel 
shows that the amplitude of uiag, Aag, decreases with increasing 
redshift width, as one would expect [equation <12H . However, A^s 
does not follow the 1 / Az behaviour predicted | Bouche et al. 200^ 
and equations <1 H & <12H . This is easily explained given that Az is 
a selection criteria upon photometric redshifts: doubling Az does 
not mean we doubled the width of the (true) redshift distribution 
since the finite uncertainty in the photometric redshifts, az, is com- 
parable to Az. 

As noted in Section l4m the bias-ratio 6]yiaii/&LRG is equal 
to a in the case of a top-hat redshift distribution N{z). However, 
given the uncertainties in the photometric redshifts, our LRG sam- 
ple is distributed arou nd the abs orber redshifts in a G aussian man- 
ner. In lBouche et al.l f2004) and lBouche et alj i2005l> . we showed 
that, in the case of a Gaussian redshift distribution, the amplitude of 
the Mg 11-LRG cross-correlation relative to that of the LRG-LRG 
auto-correlation, a, is overestimated by 25 ± 10 per cent. This cor- 
rection factor was determined using (i) numerical integrations and 
(ii) mock catalogues (from the GIF2 collaboration. lGao et alj200^ 
made of galaxies that had a redshift uncertainty equal to the slice 
width, Az, as in the case of our LRG sample. Thus, the bias-ratio 
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Figure 7. The solid line shows the bias b(Mi^) as a function of the halo- 
mass Afh (Mq ) using the formalism of Mo & White 1 2002 ). The thick sohd 
line shows the linear approximation over the mass range log A/jj = 11.8- 
13.2 used to estimate that our Mgll absorbers have haloes with mass 
(log A4(Mo)) = 11.941q [equation |^] from the bias ratio [equa- 
tion 1221 1. The dotted line shows the linear approximation over the mass 
range log A/ijCMq) = 10.7-13.5 used in Section l4!6l The inset zooms in on 
the region of interest. 



inferred from the measured amplitude ratio equation <21H is 



"Mgll 

&LRGs 



= 0.65±0.08(stat)±0.05(sys), 



(22) 



ms the un c 

comparison, in Bou che et alj 12004'), we used 212 Mgll absorbers 
and ~ 20, 000 LRGs selected in SDSS/DRl and found that the am- 
phtude ratio was 0.67 ± 0.09. 

This bias measurement is entire ly consistent w ith the results of 
iBere eron & Boisse 1 1991) and Stei del et ai]<1994 who found that 
Mgll absorbers with W,^'^''^^ ^ 0.3 A are on average associated 
with late-type ~ 0.7Lg galaxies, since the expected amplitude ra- 
tio be tween early and late type galaxies is ~ 0.70 (see Bouche et al. 
|2004) . In addition, the auto-correlation of late-type galaxies has a 
shallower power-law slope than that of re d early-type galaxies (e.g. 
iMadswick & et al.j2 003: Collister & Lahav 2005), a behaviour re- 
covered in our cross-correlation Wag measurement in Fig.|3l 



4.5 Halo-masses of Mg ii absorbers 

As already mentioned, in the context of hierarchical galaxy forma- 
tion, equations (S) & |9j imply that equation <22t can be used to 
infer the halo-mass of the Mg II absorbers provided that the mass 
of LRGs is known, which is the case here, as described in Section 

However, the transformation of equation <22t int o a halo-mass 
is no t entirely straightforward. Indeed, as noted in iBouche et alJ 
i2005i) . equations jsj & |9) refer to the mean bias, b, averaged over 
the mass distributions fe(> A/) = p{M') b{M') d\og M' , 
where p{M) = ^ is the appropriate mass distribution [nor- 



malized such that J p{M) dlog A/ = 1] and b{M) is the bias of 
haloes of a given mass M, shown in Fig.0 In the case of the LRGs, 
p(A/) is given by the halo mass function n{M). Unfortunately, the 
distribution p{M) = dn/dz/dlog A/ for the Mgll absorbers is 

unknown. 

As discussed in lBouche et aljiioOSh . this issue may be allevi- 
ated by expanding b{M) to first order over a restricted mass range, 
i.e., b{A'I) = bo + bi log M. The mean bias b is then given by 

/•oo 

6(>Af„,in) = / p{M')b{M') dlogM' 



bo + bi 



p{M') logA/'dlogAT' 



= bo+bi{logM) , 



(23) 



where (log Af) is the first moment of the distribution p{M). For 
example, the mean bias of the Mgll absorbers is 6 evaluated at 
the mean halo-mass, i.e. ^jyigu = &((log Afh)), where <> is the 
average using the appropriate mass distribution p{M). Since the 
coefficients bo & 6i can be obtained from the bias function b{M), 
the distribution p{M) need not to be known. 

The thick dashed line in Fig. shows the linear fit to 6(Af ) 
over ~ 1.4 dex. Provided that the mass range of interest is small 
(< 1 dex), this approximation is valid. Using this linear fit to b{M), 
our bias-ratio measurement equation <22t implies that the Mg 11 ab- 
sorbers reside in haloes of mean mass 



(logA/h(MQ)> = 11.94 ± 0.31 (stat) 



(sys), 



regardless of the actual p{M) distributions. An additional source 
of systematic error is the mass of LRGs, which has been assumed 
to be ~ 3 X 10^'' Mq. Given that the uncertainty in the correlation 
length ro,gg is of the order of ~ 1 Mpc at most (Section |4.2> . we 
find that the uncertainty in the LRG halo-mass is Iq ^s dex. This 
implies that the additional systematic uncertainty in the absorber 



halo-mass is ' 



dex. 



By adding this additional systematic uncertainty in quadra- 
ture, we find that the mean halo-mass of Mg 11 absorbers is (with 
1-(T uncertainties): 



(logAfh(M0)> = 11.94 ±0.31(stat)l!3;^^ (sys). 



(24) 



Adding the statistical and systematic errors in quadrature gives 
(logA/h(M0)> = 11.94l°:^g. 

Is this mass scale ~ 10^^ Mq consistent with the incidence 
probability of Mg 11 absorbers? In other words, are there too many 
(or too few) such haloes in a ACDM universe? The incidence prob- 
ability of Mg 11 absorbers, dN /dz, is given by the volume number- 
density of haloes, n, times the co-moving cross-section, tJco: 



dN 
d7 



dr 
dl 



(25) 



where the cross-section CTco ~ ^fipiiys(l + -^^abs) for 
i?phys the radius of the cr oss-secti o n in physical units. 
Since i?phys is ~40ft~^kpc jSteidell \l99^ and the num- 
ber density of h aloes of mas s M = 10^^ M© IS n = 
lO-^/i"^ Mpc-3 <Mo&Whitj|2002 h. we find that dN/dz ~ 
0.3 (n/10^^) (-Rphys/40 kpc)^. This is cl ose to the observed 
value for Mg II a bsorbers with W^^'^'^^'^ > 1 A jNestor et alj2005t 
IProchter et all2005i^ . We can therefore conclude that there are nei- 
ther too few nor too many Afh = 10^^ Mq haloes to account for 
the observed incidence probability. 
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Table 3. The bias-ratio and the inferred halo-mass of Mg II absorbers, with 
1-(T statistical uncertainties, as a function of W^^"^^^, binned in two differ- 
ent ways. The first splits the sample of nabstot Mg II absorbers into 4 bins of 
equal size except for the last bin which is double the size to provide enough 
statistics. The second uses 3 bins of equal size. Regardless of the binning, 
the bias-ratio decreases with increasing W^'^'^^^'. 





a = «)ag/Wgg 


^Mgll/^LRG 


logAf]y[gIl{MQ) 


0.3-5.5 


0.808 ±0.096 


0.65 ± 0.08 


11.94 ±0.31 


0.30-1.15 
1.15-2.00 
2.00-2.85 
2.85-5.50 


0.99 ±0.11 
0.81 ±0.11 
0.54 ±0.17 
0.66 ± 0.24 


0.79 ± 0.09 
0.65 ± 0.08 
0.43 ±0.13 
0.53 ±0.19 


12.49 ±0.34 
11.93 ±0.33 

11.11+0-52 

11.47+0.74 


0.725-1.575 
1.575-2.415 
2.415^.115 


1.13 ±0.12 
0.70 ±0.13 
0.61 ±0.14 


0.89 ± 0.09 
0.56 ±0.11 
0.49 ±0.11 


12.88 ±0.37 
11.59 ±0.41 
11.30 ±0.43 



4.6 Are Wr and Mh anti-correlated? 

We searched for possible correlations between the halo-masses of 
our 1806 Mgll absorbers and various parameters of our survey 
such as the absorber redshift Zabs, the quasar redshift Zqso, Mgll 
equivalent width, signal-to-noise of the SDSS spectra, W,^'^''^^, the 
quasar magnitude and the quasar colours. We did not find any sig- 
nificant correlation with these param eters except for the equiva- 
lent width, Wr-^2796 Qj^j sample jBouche et all 2004') already 
showed tentative evidence that the cross-correlation had a slightly 
smaller ampfitude for higher equivalent width absorbers. Here we 
significantly strengthen that evidence. 

Figure Isjleft) shows the bias-ratio &]yigjj/&LRG for sub- 
samples in 4 equivalent width bins of equal width except for the 
last bin which is twice as large in order to have enough statistics. 
This figure shows that the bias-ratio, and therefore the halo-mass, 
is anti-correlated with W^'^'^^^. We also divided our sample into 
3 sub-samples of 602 absorbers. Table|3|lists the halo-masses for 



each W, 



' bin and shows that the anti-correlation is robust under 



different binning. The bias-ratio in this figure is computed from the 
ratio of the Mg II-LRG cross-correlation to the LRG-LRG auto- 
correlation measured on large scales (1-13 h^^ Mpc) as in Sec- 
tion 4.4. In other words, the anti-correlation shown in Fig. 8 has 
nothing to do with 'Mgll absorbers in the haloes of LRGs', since 
the haloes of LRGs are much smaller than the scales probed by our 
clustering. 

Figure |8[right) shows the halo-mass inferred for each of the 
4 bins used in the left plot. In converting bias-ratio to halo-mass, 
we used the same method outlined in Section R31 However, since 
the bias-ratio in Fig. |8jleft) corresponds to a larger mass range 
than in Section 14.51 we adopted the dotted line shown in Fig. |7| 
as a linear approximation to the bias function [equation )23H . This 
approximation is naturally poorer than for our mass estimate for 
the entire sample [equation <24H . Of course, equation <23> can 
be expanded to the second order, but it would include the second 
moment ((log A/)^) of the distribution j3( A/), which is unknown. 
Quantitatively, we find that absorbers with W^'^'^^'^^ ?J 2 A have 
haloes with (logAfh(MQ)) = ILSlg*, while Mgll absorbers 
with 0.3 < M/,^2796 < 1 jjg^g (logAfh(MQ)) = 12.5lo;3, 
i.e. a mass difference of > 1.0 dex. 

The observed Afh-W^^^^^® anti-correlation in Fig.|8|has im- 
portant implications. If Mg II clouds were virialized entities in the 
gaseous haloes of the host-galaxies, one would expect that the ve- 



locity spread of Mg II systems would be proportional to the mass of 
the host galaxies. At equivalent widths Wj^2796 ^ q g jj 
transitions of the sub-components are saturated and so 14^/^2796 
strongly correlates with the velocity spread (Av) of the absorber 
(e.g. Ellison 2006). This velocity range for strong Mgll systems 
covers the range 50^00 kms~^. If the individual clouds are viri- 
alized in the halo of the host-galaxies, Av ought to be related 
to the velocity dispersion of the gaseous halo and to the mass of 
the host galaxy. Therefore, the M-t^-W^'^'^^^ anti-correlation shows 
that most strong Mg II absorbers are not virialized in the gaseous 
haloes of their host-galaxies. 

In fact, it is often tacitly assumed in the literature that Afh 
and Av (as traced by W^'^'^^^) should be positively correlated and 
so the M-ti-W^'^'^^^ anti-correlation in Fig. |8| may seem surpris- 
ing at first. In the next Section, we show that this auto-correlation 
is compatible with numerous past results in the literature. More- 
over, those past results actually seem to require a My.-W^'^''^^ 
anti-correlatio n, independent of our new results. 

Recently. IProchter et alj j2005li used the redshift evolution of 
AN /dz and simple cross-section arguments to infer rough esti- 
mates of the halo-mass of Mgll absorbers. Their conclusions are 
similar to ours: absorbers with large equivalent width, W^'^'^'^^ > 
1.8 A, were found to be more likely associated with Afh < Af* 
haloes, not with massive AL ~ 10^^ Mq haloes. 



5 DO PAST RESULTS REQUIRE A Mh-W^ 
ANTI-CORRELATION? 

In this Section, we study the incidence probability of absorbers 
in more detail. Specifically, we use two arguments to show that 
a M\^-W^^''^^ anti-correlation is already required by past results, 
i.e. regardless of our cross-correlation results. In Section |5!T1 we 
show that the cross-section-luminosity relationship combined with 
the cross-section-equivalent width relationship require an anti- 
correlation between luminosity and equivalent width. In Section 
15.21 we show that the cross-section-luminosity relationship, the ob- 
served incident probability (A? N / Az / AW^) and any plausible host- 
galaxy luminosity function (LF) together require an anti-correlation 
between luminosity and equivalent width. We also show that the 
slope of the expected anti-correlation is completely consistent with 
our new result in Fig.|8| 

The starting point for the discussions below is the incidence 
equation, equation J25> . but in its integral form: 



W2 



A^N 
AWAz 



{W,)AWr 



A^N 



Wi 
Ma 

Ml d log MAz 

A log A/ 



A'^N Ar 
^j-^(Af)a.o(Af)-, 



(26) 



where d^A''/dlogAf/dV is the number of haloes (or galaxies) per 
unit volume per unit mass (or luminosity). In addition, a cover- 
ing factor Cf should be included in the right-hand side of this 
equation since o bservations imply that it may be less than unity: 
iTripD & Bowei] |2005 ) studied close galaxy-QSO pairs and found 
that ~50 per cent of galaxies do not produce any Mg II absorption 
down to 0.3 A, indicating that the mean covering factor for our sam- 
ple is Cf ~ 0.5. We note that Cf is likely to depend on the size of 
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Figure 8. Left: The bias-ratio 6(Ma)/fe(Mg) as a function of rest-frame equivalent width 



using 3 equally spaced bins plus a larger bin. The bias- 



ratio is computed from the ratio of the Mgll-LRG cross-correlation to the LRG-LRG auto-correlation measured on scales 1-13 Mpc as in Section 4.4. 
The bias-ratio is anti-correlated with W^^""^^^. This anti-correlation implies that higher equivalent width absorbers are less massive. Right: The halo-mass, 
Mij(Mq), of Mg II absorbers as a function of W^^'^^^ in bins corresponding to those on the left plot. Absorbers with W^^"^^^ > 1.5 A have a halo-mass 
1 dex smaller than absorbers with 0.3 < VFi'^rge < i,o A. 



the host galaxy and therefore its mass, a comphcation that we ig- 
nore here. 

For the past two decades, the main problem in solving equa- 
tion <26> for the cross-section CTco was that one had to perform the 
integrals of the left-hand side over a range of W/'^^^® and the right- 
hand side over a range of mass M (or luminosity L), ignoring any 
possible - and probably quite strong - dependence of W^'^"^^^ on 
L or M. We will show in the next two sub-sections that the A/h- 



anti-correlation is required (independently of our cross- 



correlation results). 

We first show (in Section 5.1) that a Mh-Vl^r^^^^^ relationship 
is required from the following simple arguments. Observationally, 



for a given bin in 



there is a maximum impact parameter, 



^phys, that defines the cross-section, a. The radius of the cross- 
section depends on the luminosity of the host galaxy, iiphys oc L'^, 
i.e. the cross-section is a function of mass or luminosity. Since 
there is a relationship between W^^'^^^ and iiphys and between 
i?phys and L, there is a relationship between W^^'^^^ and Mh (or 
L). Section 5.1 shows that the two physical parameters should be 
anti-correlated. Section 5.2 arrives at the same conclusion based on 
different arguments which make one assumption, namely that the 
relationship between log Mh and W^^'^^'' is linear. 

5.1 The cross-section distributions require a Mh-W^^^^^ 
anti-correlation 

Here, we investigate the existence of a luminosity (or mass)- 
^/^A2796 relationship and whether the two parameters should be 
correlated or anti-correlated given the observed distribution of 
Mg II absorbers in the cross-section-equivalent width plane and the 
dependence of the cross-section on luminosity ( Steidel 1995). 

The radius of the cross-section, i?phvs, decreases steepl y with 
increasing W^^'^''^'^: Lanzetta & Bowen ( 1990) and Steidel Il99^ 
(their figure 3) showed that the host-galaxies of absorbers with 
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Figure 9. The Mgll equivalent width, W^'^'^'^^ , as a function of impact 
parameter, p, for the host-galaxies of absorbers with W^2796 > q 3 ^ 
The solid squares show the impact parameters for host-galaxies selected 
from a sample of our SDSS Mg II absorbers closer than 100/i~^ kpc (co- 
moving) to the QSO lines-of-sight and with jzphot ~ -^absl < 0. 05. The 
open squares, op en triangles a nd s olid triangle show the samp les of lSteidej 
>1995ll) . lKacDrzaketalJ (2006) and B ergeron & Boissj 11991B respectively. 
All impact parameters are shown in our cosmology. The solid line shows 
the fiducial model for the maximum impact parameter Rphys used in Sec- 
tion lS.lK see text) to represent the cross-section-equivalent width relation- 

ship, CTcoiW,^^^^<^). 
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^^A2796 > 2 A do not exist at large impact parameters. Figure |9| 
reproduces these results. The data points show the impact parame- 
ters p (in physical units) of the Mg II host galaxies from the samples 
oflgeraeron & Boisse 1 1991) (filled triangle), Steidel ( 1995) (open 
squares) and lKacprzak et aU ^2006) (open triangles) for absorbers 
with W^^^"^^^ > 0.3 A. The solid squares represent the 44 host- 
galaxies that we have identified in the SDSS/DR3 imaging data 
with l^abs — ^phot] < 0.05 and which were within 100 kpc (co- 
moving) of the QSO line-of-sight. This sample is truncated at an 
impact parameter p < 20 kpc (physical; ~3") due to the large 
point-spread function of SDSS. It is clear that the data points do 
not fill the upper right side of the plot. More precisely, at increasing 
equivalent width, the maximum impact parameter ~ which defines 
the radius -Rphys of the cross-section - decreases. This dependence 
can be roughly parametrized as Wr(-Rphys) — —3 log[iiphys] + K 
as represented by the solid line in Fig.|9| 

In addition, a relationship between between -Rphys and lumi- 
nosity L also exists: ISteidell jl995h found that -Rphys slightly in- 
creases with luminosity: 7iphys(I') oc {Lk / Ll^)^ , with /3 ~ 0.2. 
Thus, the simple observational results that there are relationships 
between Wr & -Rphys and between -Rphys & L require a L-W^^''^^ 
relationship. 

We emphasize that this L-W^'^'^^^ relationship does not mean 
that there is a direct correspondence between mass and equivalent 
width on a galaxy-by-galaxy basis. This relationship simply means 
that on average a sample with a well-defined W^'^'^^'' (e.g. a bin of 
^^A2796^ will have a well-defined and predictable mean mass (or 
luminosity). More specifically, for a sample of absorbers selected 
in a bin of W^^'^^^, the host galaxy will have an impact parameter 
less than -Rphys, and a mean mass A/h- One should not interpret the 



distribution in 



as being entirely due to the distribution in 



halo mass (or in galaxy luminosity). 

Finally, if one combines the observed cross-section- 
equivalent width relation (Fig. |9) with this cross-section- 
luminosity relation one finds that W^^^" (log Lk) oc 
-3/3 log Lk ^ -0.6 log Lk, i.e. log Lk oc ~ -l.GWr- 

In other words, the cross-section as a function of luminosity and as 
a function of equivalent width together require that Wr (therefore 
Aw) and luminosity (or mass) be anti-correlated. 
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Figure 10. The sym bols represent the absorber K-hwA luminosity function 
rLFiof lSteideiJl99 .5V The error bars are uniform and represent the average 
uncertainties. The solid line represents a Schechter function w ith a = —1, 
as observed for the field -ft'-band LF iKochanek et alJIlOOll) . The dashed 
line represents a Schechter function with a = —0.5. Both lines are reason- 
able descriptions of the -ftT-band LF. Both curves were scaled to match the 
data points. 



An independent way to show the existence of a L-Wr anti- 
correlation makes use of equation <27> in terms of luminosity: 



dzdWr 



dlogL _d^/V_ 

dWr 



dr 
dz 



(28) 



dVdlogL 

where d^-ZV/dV^/dlog-L is simply the LF in logarithmic 
form, 0(L)L In(lO). The differential incidence probability, 
d^N/dz/dWr, is now well co nstrained observationally thanks to 
the la rge SDSS database (e.g. iNestor et al]l2005l : IProchter et"ai] 
l2005t) . 

For luminosities below L*, (f){L)L oc L^'^°', where a is the 
faint-end power-law slope of the LF. Equation <28t then implies 
that 



5.2 The incidence probability requires a Mh-W^ 
anti-correlation 

Since we have established the existence of a M^^-W^^ 
ship, the right-hand side of equation <26> can be rewritten as 

rW2 



' relation- 



dW, 



dlogM d^iV 



Wi 



dWr dlogMdV 



dr 

dz 



By differentiation of both sides of equation <26> . it follows that ' 



d^N 
dWAz 



dlogM d^iV 



dW, dlogMdl/ 



dr 
dl 



(27) 



where ',°,f/^ is the Jacobian of the transformation between Wr 
and M. We will assume that the relationship between Wi and AI is 



linear, i.e. 



d log M 
dWr 



IS a constant. 



3 If F(x) = /(a)da and G{x) = J"" g{a)da. If F{x) = G{x) for 
all X then F = G and their derivatives are also equal, /(a) = g(a). 



d^N 



r(l + a + 2/3) 



^ ^1 T - " (29) 
dzd log L 

using (Tco(-t') oc L^^ from lSteideilh995h . In addition. iNestor et all 



mO$) showed that d^N/dz/dW, is an exponential, i.e. 

exp(-W./H/.), 



dzdW, 



(30) 



where Wt is the exponential scale length. 

Quantitatively, from equations <29> and <30> , the expected 
slope of the L-Wr anti-correlation is 



d log L 

dWr 



-log,o(e)/W^* 
l+a + 2l3 ' 
-1.8 fora = -1.0, or 
-0.8 for a = -0.5, 



(31) 



since INestor etalJ j2005l) showed tha t W* ~ 0.6 i n the redshift 
range (z) ~ 0.5 of our survey, and lSteidell il995|) showed that 
(3 = 0.2. 

Therefore, equations <29> and <30t imply that Wr and log L 
(or log Ad, using a constant M / L) must be anti-correlated as long 
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as 1 + a + 2/? > 0, i.e. a > -1.4 for (3 = 0.2. Most plausible 
LFs easily satisfy this condition. 

The S-band LF of Mgll absorbers was constrained by 
ISteidel e t al. 1 1995) and is known to be bell-shaped, with a shal- 
lower than -1 (figure 4 of Steidel 1995). The if-band LF of Mgll 
is likely to be closer to a mass function. Figure [Tol reproduces the 
absorber iC'-band luminosity function of Steidel ( 1995). The er- 
ror bars are uniform and represent the average uncertainties. The 
solid line represents a Schechter function with a — —1, as ob- 
served for the z = ii'-band field -galaxy LF iKochanek et alJ 
l200lHBell et all2003h . lSteidej Jl995l ) argued that the host-galaxy 
if-band luminosity function is similar to the /("-band field-galaxy 
LF. The dashed line represents a Schechter function with a — —0.5 
and is an equally good match to the data. 

One can go a step further and 'predict' the Mh-W^^'^^^ re- 
lationship using the mass-to-light ratio. Since the faint-end slope 
of the mass function of haloes (e.g. Mo & White 2002) is steeper 
than that of galaxies, t he halo mass-to -light ratio, Mh/L, decreases 
with increasing L (see lShankar et al.t2006. and references therein). 
Parametrizing Mh/L as L'' (with q < 0), one finds that the ex- 
pected slope of the M^-W^^'^^^ relationship is: 

dlogM .dlogL 

= (1 + g) 



Vc/c). However, as demonstrated in Sections 5.1 and 5.2, one 
would not reproduce the observed relation between i?phya and 



. Thus, both our cross-correlation results and the arguments 



dW,. 



+ 



dWr 

1 + a + 2/3 



(32) 



Fig. II If left) plots the expected relationship between Afh and 
W,^^'^^^ for g = 0, -0.5 & -0.75 and a = -1.0 & -0.5, with all 
curves normalized at log Mh (Mq ) = 12. The data points are as in 
Fig.|8jright) and agree well with the expected relationships (lines). 

Our conclusion from this ex ercise is tha t any plausible LF, 
combined with d^N/dz/dW. (e.glNestor et ali2005 ) and with the 
observed i?phys(iif ) relation (iStei deil995h . together require a L- 
W^'^'^^^ (or Mh-W^'^'^'^^) anti-correlation. 

Since a M^-W^^'''^^ anti-correlation is required from past re- 
sults, we can perform a simple consistency check as follows. We 
use the series of relationships shown in Fig. ll If left) and the halo- 
mass function, n(Mh) = dNIdVId log M. from Imo & Whitd 
i2002h in equation <27> to predict the radius i?phys of the 
cross-section as a function of equivalent width, R-physiW^^"^^^), 
that is required by the incidence distribution d^A^/d2:/dWr of 
iNestor et alj J2005h . The result is shown in Fig. ll IT right) for each 
of the relationships on the left plot. The data points are as in Fig.|9| 

Thus, we have shown with two different arguments (in Sec- 
tions 5.1 c& 5.2) that a positive Mh-W^'^''''"^ correlation is not con- 
sistent with the impact parameter distribution {p-W^^"^^^ ) of cur- 
rent Mg II samples. This conclusion is completely independent of 
our clustering results presented in Section|4] One can turn these ar- 
guments around and state that the anti-correlation between W^'^'''^^ 
and maximum impact parameter (Fig.|9} is a natural outcome of the 
Mh-W^/'^^^'^ anti-correlation in Fig.lH 



6 PHYSICAL INTERPRETATION 

The Adh-W^^"^^^ anti-correlation has a direct and important con- 
sequence: As described in the introduction, if the individual sub- 
components or clouds of strong Mgll systems were virialized 
within the host-galaxy halo, the velocity spread Ad (as traced by 
W^,-^2^^^) would be a measure of the galaxy velocity dispersion and 



would be correlated with Mh via W,. 



7A2796 = Av/c (or 



with the circular velocity of the galaxy, 14, via ^''/A2796 



presented in Section 5 show that the Mg II clouds are not virialized 
in the gaseous haloes of the host-galaxies, in the sense that the ve- 
locity spread Av (as traced by W^^'^^^) of the sub-components is 
not proportional to halo mass. 

We now turn towards a physical interpretation of our results. 
In particular, can we put constraints on the competing models for 
Mg II absorbers? 



6.1 Basic Mg ii cloud properties 

In a few cases the physical properties (sizes, densities etc.) of the 
individual components of Mgll systems with W^^^^^'^ > 0.3 A 
have b een constra i ned f rom photo-ionization mod elling. For in- 
stance, IBin^ni] j^OO?) and Masiero et alj j2005h found that in- 
dividual Mgll velocity components (or clouds) of systems with 
i[^A2796 ^ Q j^jjyg densities of nn = 0.001-0.1 cm~^ 
and, therefore, typical siz es of s = 10-1000 pc, close to the di- 
rect size measurements of lRauch et alj i2002h using lensed QSOs. 
For a warm (T ~ 10* K) cloud of size s = 100 pc and nn ~ 
0.01 cm"'', the sound crossing time is ~ lO'^yr, and the cooling 
time-scale is kT/[nA{T)] ~ lO'' yr since the cooling function 
A(r) ~ 10"^^ergcm^ s"^ at T = 10** K jSutherland & Dopitl 
Il993l) . The same ionization studies indicate that the associated Mg I 
components are denser (nn >10cm~^) but very small, only 0.01- 
0.1 pc. This supports the view in which the clouds would condense 
out of the hot gas. 

For such warm clouds embedded in a hot (T ~ lO'^ K) ionized 
medium, the evaporation time scale is > lO'^ yr <Cowie & McKed 
I19771 Im cKee & Cowie '1977^. They are stable against gravity be- 
cause their leans lengths are much larger than the cloud sizes, but 
they would be destroyed by sho cks in s/Vshock^ 1 Myr for a shock 
velocity fshock of lOOkms"^ iKIein et alll994 . 



6.2 HVCs and Mg ii clouds 

The properties of the Galactic high-velocity clouds (HVCs) - 
their H I column densities, sizes, densities e tc. - are akin to those 
of Lyman limit and Mg II absorbers (e.g. [ Charlto n et all l200Ct 
Zwaan & Briggs 2000). For instance, iRichter 6^1712003) stud- 
ied the Galactic gaseous halo in high-resolution spectra towards 
one line-of-sight (PKS 1448—232) and detected low column den- 
sity, high-velocity Call (iVcaii ~ 10"~^^cm"^) and Nal D 
(A''Nai ~ 10^^ cm"^) absorption at VhSR — — 150kms"^ with 
Doppler widths of 4-8 k ms~^. This is very sim ilar to typical val- 
ues for Mgll absorbers I Churchill et al J l2003h . Their follow-up 
VLA 21 -cm studies unveiled several high- velocity Hi clumps at 
peak column densities A^'hi ~ lO'^^ cm~^ near the QSO line-of- 
sight. The sight-line seems to have intercepted the outer part of 
one of the clumps. Assuming a distance of 4-12 h^^ kpc, the den- 
sity of this clu mp is nn — 0.1 c m~'', similar to that inferred for 
Mg II clouds bv lDing et alJi200'3l) . Because 8 of 13 sight-lines have 
similar high-velocity Call clouds. iRichter e t al. (2005) concluded 
that the covering factor of such low column density gas must be 
large. Based on this they argued that if such clumps are typical for 
haloes of quiescent spiral galaxies, such low column density clouds 
should contribute significantly to the population of Mg II absorbers 
and Lyman limit systems. In a starburst environment these clouds 
may be either more numerous and/or have larger column densities. 
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Figure 11. Left: The black squares show the halo-mass (Mh) as a function of equivalent width ( W^^'^^^) as in Fig.|8jright). The grey squares show a different 
binning. The various Unes represent the expected Mi^-W^^"^^^ relationship, normalized at logMji(MQ) = 12, for various LFs (a = —1.0, —0.5) and 
mass-to-light ratios (q = 0, —0.5, —0.75) [see text and equation I32i 1. The results from our cross-correlation analysis follow the expected anti-coiTelation. 
The error bars represent the l-cr statistical uncertainties and the size of the equivalent width bin used. Right: W^'^'^^^ as a function of impact parameter, p, 
for the Mg II host-galaxies. The data points are as in Fig.|9] The lines are as in the left panel and represent the radii of the cross-section J?phys (i.e. maximum 
impact parameter allowed). The numbers along the dot-dash line indicate the logarithm of the halo mass corresponding to the W^^"^^^ from the left panel. 



Towards the Andromeda galaxy fM3 11. iThilker et alJ 
reported a population of faint Hi clouds detected directly 
with the G reen Bank Telescope (GBT). Follow-up studies by 
IWestmeieTet alj i2005l) using the Westerbork telescope clearly 
show that the clouds have typical sizes of ~ 1 kpc, Hi column 
densities of 10^^"^" cm~^, densities nu ~ 10^^ cm^'', Hi line 
widths of ^ 20kms~^ and velocity gradients of several kms"^. 
The Hi mass of these clouds is then '--^ 3 x 10^ M0. Thus, the 
properties of HVCs are very close to those inferred for Mg II clouds 
from photo-ioniz ation studies disc ussed in Section lSm Csee also the 
recent analysis o f iFox et alfcOO^) . 

6.3 Finding the right model 

There are three main hypotheses that could describe Mgll ab- 
sorbers. The first assumes that QSO sight-lines probe galactic 
gaseous disks and, therefore, the interstellar medium of the host- 
galaxies. The second assumes that each QSO sight-line probes 
gaseous haloes of galaxies where the gas originates from the fil- 
aments, falls into the galactic potential well and eventually cools 
onto a rotationally supported disk. The third hypothesis places the 
Mg II absorbing gas in the same gaseous halo, but the gas comes 
from outflows driven by the large number of supernovae in star- 
bursts. This last hypothesis has long been used to model the HVCs 
in the conte xt of a 'galactic fountain' scenario in a m ulti-phase 
medium (e.g. lShapiro & FielJl97ilWolfire et alll995h . 

The disk model predicts a strong dependence of W^^^^^ 
on inclination of the disk with respect to the line-of-sight. 
[iCacDrzak et al. ( 2005, 2006} did not find any such correlation. Fur- 
thermore, ISteidel et alj <2002t) found that, in a small sample of 5 
Mg II absorbers, an extension of the host-galaxy rotation curve did 
not reproduce the absorption-line kinematics. A simple disk model 
therefore seems unlikely. 



The evolution of the Mg II incidence probability is observed 
to strongly mimic the evolution of the star formation rate density 
(SFRD) with redshift, as pointed out by Prochter et al. ( 2005) and 
also Misawa et alj 1200 5). Indeed, removing the cosmological de- 
pendence (i.e. dr/dz) from the incidence probability, dA''/d2:, in 
equation <27t leads to the line-density of absorbers per co-moving 
Mpc, dN /dX = n{M) CTco- This line-density shows a decline at 
z < I, just like the SFRD, and both dN/dX and the SFRD are 
approximately constant from z ~ 3 to 2: = 1. The decline of 
dA*'/dX can only be explained by a decrease in the cross-section, 
(Jco, with time since the number density of haloes, n {M), increases 
with time for all masses. From these arguments, IProchter et alJ 
(2005) concluded that "the processes responsible for strong Mgll 
absorbers are turning off at 2; ~ 1". This behaviour is natural within 
a starburst scenario where Mg II absorbers are directly related to 
supemovae-driven winds. 

R egarding the in-fall hypothesis. Mo & Miralda-Escudd 
Jl99fl) produced a model for Lyman limit systems, Mg II and C IV 
absorbers in a CDM cosmology. This model consists of a two phase 
medium for the gaseous halo. In the context of galaxy formation in 
dark matter haloes (e.g. White & Rees 1978; White & Frenk 199 1*), 
the initial gas in a galaxy (either at the time of collapse or through 
accretion) is shocked at around the virial radius. A halo of hot gas 
(T ~ lO*" K) is formed which has a temperature close to the halo 
virial temperature since it is in an approximate hydrostatic equilib- 
rium with the dark matter. Due to cooling from thermal instabilities 
and inhomogeneities, a warm phase (T ~ 10'' K) will form, com- 
prising photo-ionized clouds confined by the pressure of the hot 
phase, while accreting onto the galactic central region. 

As a consequence, in this inflow picture the photo-ionized 
clouds will be appro ximate ly virialized in the gaseous halo. 
|Mo & Miralda-Escud^ il996h argued that the terminal velocity of 
these clouds will be close to the virial velocity for haloes with 
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Vc < 250 kms"^ In this scenario, the velocity spread of the gas 
All (as measured by W^^^^^) would be related to the virial veloc- 
ity and therefore its halo-mass. Thus, we view the observed Afh- 
p[/A2796 anti-correlation as strong evidence against this inflow pic- 
ture. Furthermore, more massive (and larger) galaxies with large 
W^^"^^^ would be seen at large impact parameters, contrary to th e 
observations of Lanzetta & Bowen ( 1990) and Steid el et alJil99^ . 

There is on e main caveat to this conclusion. The in-fall model 
of IMo & Miralda-Esjaide 1 ,1996) might not be directly relevant 
to our absorbers with W^^""^^^ ~ 1-4 A which are expected to 
be optically thick, i.e. saturated. The in-fall model assumes the 
clouds are optically thin, i.e. W,^^''^'^ < 0.3 A, still on the lin- 
ear part of the curve of growth. Most such absorbers are pre- 
dicted to be in massive galaxies, with a median Vc ~ 220kms~^ 
[logMh(M0) ~ 12.4] and most systems in the range Vc = 150- 
300 kms"^ [logMh(M0) ~ 11.9-12.8]. This mass range, pre- 
dicted for 0.3 a absorbers, agrees with our mass esti- 
mates for the lowest W^^"^^^ systems in our sample (0.3-1.0 A). 
The inflow model could therefore describe Mgll absorbers with 

^/A2796 < -^g gj^^, jjpjg jjjgj jjjg ^;-A2796 (Jistj-ibution 

function has different power-law indexes below and above 0.3 A 
(Nestor et al. 2005), possibly reflecting two different physical ori- 
gins. 

Another possible caveat is the recent realisation that, for low 
mass galaxies, a significant fraction of the gas accretion occurs 
without the initial shock heating because the cooling time is very 
short, i.e. the gas remains warm (T ~ 10* K) as it moves from 
the intergalactic medium to the central parts of galaxies. For this 
reason, it is often referred to as the 'cold mode' of gas accretion 
(e.g. Birnboim & Dekel 2003; Keres et al. 2005). In this context, 
the clouds will not be virialized, but the terminal velocity should 
still be related to the ha lo -mass. 

^__|Mallej^^Bullgc^ J2004|) extended the arguments of 
IMo & Miralda-Escudd i 1996ft to include multi-phase cooling 
since the hot gas in galactic haloes is thermally unstable (see 
also discussion in Section l6m . This multi-phase cooling scenario 
allows the survival of a hot gas core and naturally gives rise to the 
formation of ~ 10'* K clouds, each of mass ^ 10*^ Me, which 
are pressure-confined by the hot halo. The main conclusion of 
this model relevant to our work is that this multi-phase treatment 
naturally explains the properties of Mgll clouds and HVCs: 
the unstable hot gas produces warm clouds of mass ^ 
which are stabilized by the pressure of the surrounding hot gas. 
In this model, the maximum distance or impact parameter for 
warm clouds (HVCs) is set by the the 'cooling radius', which is 
'-^ 100 kpc for a Milky Way size halo. This model differs from 
the one of Sternbere et alj i2002h in that the self-gravity plays a 
negligible role and no dark matter is required. 

IStemberg et aL (2002) modelled the HVCs as multi-phase 
(warm-ionized and cold-neutral) pressure-confined cloud, each be- 
longing to a 'mini-halo' made of dark matter. In this model, the 
dark matter in the mini-haloes provide the extra pressure in the 
core of the clouds required to allow for the formation of a cold 
medium. The warm phase of the HVCs is confined by the am- 
bient pressure provided by the hot (T ~ 10® K) ionized gas 
of the halo. This hot gas is in hydro-static equilibrium with the 
dark-matter of the host galaxy and determines the pressure profile 
of a galaxy. Given the observed velocity widths (10-20 kms~*), 
IStemberg etai]<2002ft determined that the pressure of this hot-gas 
must be P/k > 50 cm"'' K to keep the warm phase of the clouds 
bound. In the Milky Way, the pressure of the hot halo is above 
50 cm"'' K out to 150-200 kpc, i.e. almost to the virial radius. 
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Figure 12. Equivalent width, , as a function of impact parameter, 

p, as in Fig. II ir rightl The data points are as in Figs.l9l& llll ri ght) . The 
solid line shows one of the model curves from Fig. II II the long dashed line. 
The numbers along the solid line indicate the logarithm of the halo mass 
in units of Mq) corresponding to the equivalent width. The dashed line 
shows the virial radius of such haloes (r2oo). with the corresponding virial 
velocity marked. The virial radius is much greater than the cross-section 
radius, i?pjiys , for absorbers with W^^''^^^ > 1.5 A. Therefore, strong Mgll 
absorbers seem not to fill the haloes of their host-galaxies. 

For galaxies with smaller halo-mass (or smaller virial temperature), 
the central pressure will be reduced and the radius at which P/k 
reaches 50 cm^'' K will be smaller than the virial radius, at which 
point the cloud will evaporate. 

The fact that no absorbers with large W^'^'^^^ are found to 
have host-galaxies at large impact-parameters (i.e. no points to the 
right of the solid line in Fig.|9} clearly tells us something about 
the characteristic scale where clouds either form or evaporate, and 
not the size of the halo of the host-galaxy. Indeed, for the halo- 
masses found in Fig. II If left) one can compare the correspond- 
ing virial radii to the p-W^^^^^^ distribution in Fig. II ll right). as 
we show in Figure fT2l The dashed line represents the z — 0.5 
virial radius, r200! for the model with a = —0.75 and g = in 
Fig. II ir ieft). The logarithm of the halo-masses (in Mq) are marked 
along the solid line. The halo circular velocity, Vc, is marked along 
the dashed line in kms^^. This figure clearly shows that r2oo is 
much larger than the cross-section radius (i?phys, solid line) for ab- 
sorbers with VK,-^2796 > ;^ 5 ^ (or halo-masses Mh < W^^'^ Mq). 
In other words, Mgll absorbers do not fill the halo of the host- 
galaxy. The radius of the cross-section must then be due to some 
other physical mecha nism. It could be either the pressure threshold 
discussed b y Sternber get alji2'002ft and/or to the cooling radius of 
iMaller & Bullock. ilOolt) . 

6.4 In the context of super-winds 

Local supernovae-driven winds are known to expel ionized gas 
into their haloes, as traced by Ha, X-rays and certain emis- 
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Il997h. However, the oresence of Nal D absoro- 



tion also indicates that the winds carry large amounts of dust 
and cold gas, M > lO^Mp, over the life-time of the star- 
burst (e.g.'Lehnert & Heckmaill99^ Martini 1999l:lHeckman et alj 
[2OOO; RuDke et al. 2005; Martin 20Q€ ). T he dus t can sometimes be 



seen in direct imaging ( Howk & Sa vage 2000) out to a few kpc. 
Furthermore, the amount of dust is correlated with the amount of 
cold g as traced by Na I D absorption 1 Heckma n et alj200d : lMartiril 
120061) . This shows that cold material is entrained and expelled out- 
side the disk by the i onized gas flow. In this context, the findings of 
[Kondo et alj j200d) . lYork et al. ' 1 2006) and Wild, Hewett & Pettini| 
t2006) are very interesting. For the first time, .Kondo et al, de- 
tected strong Nal DAA5891/5897 doublet absorption in two strong 
Mgll absorbers (also DLAs) at ^abs = 1.062 & 1.18 towards 
the gravitationally lensed QSO APM 08279+5255. York eTai] 
constrained the dust-reddening from the slope of stacked 
QSO spectra with foreground Mg II ab sorbers selected similarly to 
ours, but at 1 < z^hs < 1.8. Finallv. iMld & HewetJ <2005l) and 
IWild et all llOO &l found that the sub-sample of strong Mg II ab- 
sorbers showing significant Ca II absorption in the low -resolution 
SP SS QSO spe ctra are even more dust-reddened. Both lYork et alj 
and lWild et alJ found that the reddening is correlated with equiv- 
alent width, a relationship very similar to that of [Heckman et al. 
<2OO0h . 

The metallicity of our Mgll systems is unknown but 
lYork et ejl used indirect arguments to estimate average metallicities 
for a DRl sample of Mg II absorbers. They found that the metallic- 
ity decreases with increasing Mgll equivalent width^. This would 
lead to a normal mass-metallicity relation when combined with our 
Mh-W/'^'^^'^ anti-correlation. 

Given these local results and that Mg 11 absorbers predomi- 
nantly trace clouds of 'cold' material, are the Mg II clouds related to 
the entrained material in outflows ^ ? If so, then W^'^'''^^' should be 
a measure of the (past) star-formation rate and mass ejection history 
for the following reasons. Firstly, is a measure of the num- 

ber of components (or clouds), as shown bv lChurchill et alJ ilOOA 
Secondly, our results show that absorbers with W^^"^^^ > 2 A 
have smaller potential wells so the host-galaxy's supernovae-driven 
winds are more likely to eject gaseous material out of the disks. 
Thus, in the context of galactic outflows, the velocity spread (and 
therefore W^^^^^), which is just a measure of the number of com- 
ponents, is related to the (past) star-formation rate and to the mass 
ejection rate. 

Entrained clouds would have to survive from the disk to where 
they are seen at i? ~ 50 kpc. The travel time to 50 kpc is about 
lOOMyr for wind speeds of 500 km s^^. From the discussion in 
Section the formation and destruction time-scales are 1-2 or- 
ders of magnitude smaller than the travel time or halo dynamical 
time (~ 0.5-1 Gyr). It seems very unlikely that such clouds could 



^ Using a very different method to estimate A^hIi iNestor et alj i2003ft 
found a correlation between W^^'^^^' and the metallicity of DLAs selected 
with the Mg II criteria of Rao & Turnshek 1 2000) and Rao et al. ( 2006). As 
noted bv lNestor et all this relationship can be understood if higher W^^'^''^^ 
absorbers correspond to regions w ith more intens e bursts of star formation. 
^ After completion of this work, iMartinI |200(I) argued that the relics of 
cool outflows of Ultraluminous Infrared Galaxies (ULIRGs) and Lumious 
Infrai'ed galaxies (LIRGs), as traced by Nal D, will create a significant 
redshift-path density for Mgll if most L > O.IL* galaxies pass through 
a luminous starburst phase. 



survive several times lOOMyr given these physical properties un- 
less they are confined by the pressure of the hot ionized me dium. 
The same conclusion was reached by ISternberg et alj i2002l) and 
Mailer & Bullock (2004) and is discussed a bove. 

Recently, iKacprzak et alJ j2005l |2006') reported that the host- 
galaxy morphological asymmetry and W^'^'^^^ are correlated. 
While these authors favour the merger interpretation for their re- 
sults, they also discuss the possibility of star-burst driven winds. 
Again, because of the smaller potential wells of higher W^^'^^^ 
absorbers, outflows driven by spatially asymmetric star-bursts 
can create asymmetries and produce the correlation seen by 
[Kacprzak et al. ( 2005, 2006) without invoking mergers. In fact, one 
sees that most Mg II host-galaxies in their sample are relatively iso- 
lated and do not show galaxy pairs or strongly disturbed morpholo- 
gies. 



6.5 Further evidences for the super-wind scenario 

Direct evidence o f outflows exists in only a few Mg II absorbers. 
lBondetalJj200lh presented high-resolution spectra of four strong 
^^A2796 y 1 8 A) Mgii absorbers at 1 < 2 < 2 and showed that 
the profiles "display a common kinematic structure, having a sharp 
drop in optical depth near the centre of the profile and strong, black- 
bottomed absorption on either side". They interpret these features 
as a signature of super-winds arising in actively star-forming galax- 
ies. In the similar case of the z — 0.74 Mgll host-galaxy towards 
QSO Q1331-I-I7 ilEUison et all2003l) . the absorber's absolute lumi- 
nosity is quite low, L ~ 0.3L*, and the authors suggest that "the 
double-peaked absorption profiles, their striking symmetries, ... and 
large Mg II equivalent width support the interpretation of superbub- 
bles ". 

iNorman et alj J 19961) searched for Mg II absorption in QSO 
spectra that lie near the starburst galaxy NGC 520, a well- 
known local starburst with super-winds. NGC 520 has an IR- 
luminosity comparable to that of M82, shows a disturbed mor- 
ph ology and stro ng filamentary Ha emission along the minor axis. 
In lNorman et all the Mgll, Mgl, and Fell absorption associated 
with its gaseous halo are clearly detected along two lines-of-sight. 
The QSOs Q0I2I-I-0338 and Q0121-I-0327 intersect the halo of 
NGC 520 at ^ 55 h^^ kpc and ^ 25 h^^ kpc, respectively. In 
both sight-lines, the Mg II equivalent widths are greater than 0.3 A: 



= 0.33 A towards Q012I-I-0338 and W,^ 



= 1.7^ 



towards Q0121-I-0327. The latter absorber with the largest W^^^^^"^ 
has the smallest impact parameter, as expected from the larger sam- 
ple in Fig.|9| 



7 SUMMARY 

From the SDSS/DR3 we selected 1806 Mgll absorbers with 
W,^^'^^^ ^ 0.3 A in the redshift range 0.4 z^bs 0.8 and 
~ 250,000 LRGs in the absorber fields selected with photometric 
redshifts. From the ratio of the Mg II-LRG cross-correlation to the 
LRG-LRG auto-correlation, we constrained the halo-mass of Mg II 
absorbers in a statistical manner, i.e. without directly identifying 
them. To summarise our main results, we have shown that: 

(i) the ratio of the cross- and auto-correlation, uiag/wgg, is a = 
0.808 ± 0.096; 

(ii) the corresponding bias-ratio between strong Mg II absorbers 
and LRGs is by^g u/^lrg ~ 0.65 ± 0.08; 
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(iii) this bias-ratio implies thiat the absorber host-galaxies have 



a halo-mass of (log Mh(M0)) = 11.94 



+0.39 



The l-a uncertainty 



includes the statistical errors in a and possible systematic errors in 
the LRG halo-mass; 

(iv) the Mgll equivalent width, W^^'^^^, is significantly anti- 
correlated with the absorber halo-mass, Mh. 

The main consequence of point (iv) is that the Mg II clouds are not 
virialized in the gaseous haloes of the host-galaxy, in the sense that 
Av, the velocity spread of the individual clouds (ranging from 50 
to ^ 400 kms~^), does not positively correlate with halo mass. 
In other words, had the clouds been in virial equilibrium with the 
halo dark matter distribution, Av (as measured by W^'^''^^) would 
correlate with the mass of the halo. However, we find the opposite: 
lower-mass haloes harbour the host-galaxies producing the largest 
number of individual Mg II clouds spread over the largest veloci- 
ties. 

Since the Ali^-W^^"^^^ anti-correlation may initially seem sur- 
prising, we performed several consistency tests on the method itself 



and tested the Mh-Wj 



anti-correlation against numerous past 



results in the literature. For the galaxy clustering method we have 
shown that: 

• The cross-correlation is solely due to the Mg II absorbers; 

• The cross-correlation amplitude decreases when the width of 
the LRG redshift distribution is increased, as expected; 

• The Wag/wgg ratio is independent of the LRG redshift distri- 
bution, as expected if one uses the same galaxies to calculate both 
10 ag and Wgg. Therefore, the relative amplitude is free of possi- 
ble systematic errors caused by contaminants (stars or interloping 
galaxies). 

In checking past results in the literature for consistency with 



the observed Mh-Wr 



anti-correlation we found that: 



• When one combines the observed luminosity- and W^^''^^^- 
dependence of the cross-section, one finds that luminosity (or mass) 
and equivalent width must be anti-correlated; 

• When one combines the observed luminosity-dependence of 
the cross-section with the observed incidence probability of strong 
Mg II absorbers, N/dz /dWr , one finds that luminosity (or mass) 
and equivalent width must be anti-correlated for all plausible lumi- 
nosity functions with faint-end power-law indexes a > —1.4. 

In other words, a positive correlation between W^^"^'^^ and halo- 
mass (or virial velocity) - expected if the clouds were virialized - 
is inconsistent with past results. This is independent of the Afh- 
Wi-^''^^ anti-correlation fr om the clustering analysis. When one 
combines the observed M\^-W^^'^'^^ relationship with the observed 
d^N/Az/AWt one predicts a cross-section radius which agrees 
very well with the distribution of absorber impact parameters ob- 
served in previous and current samples. 

We interpret our results as a strong indication that a large 
proportion of strong Mgll absorbers, particularly those with 
^^A2796 > 2 A, arise in gal actic outflows. A conclusion similar was 
reached bv lSongail3 i2006h for z ~ 2 CiV absorbers. Supemovae- 
driven winds drive hot gas out of the galaxy plane, forming cavi- 
ties if numerous supernova explosions occur in a small region on 
a similar time-scale. The hot gas breaks out of the disk, forming 
a hot gaseous corona (traced by high-ions such as C IV and O VI) 
and may entrain cold, dusty material. Part of this gas cools, forming 
pressure-confined clouds traced by Mg II which may be the analogs 
of the HVCs. 



This interpretation allows us to make several predictions. 
Firstly, it provides a natural explanation as to why absorbers with 
i[^^A2796 ~> 2 A are not observed at large impact parameters. These 
systems should be physically closer to the host galaxy and also 
closer in time to more recent and more active star formation. Thus, 
they ought to have star formation rates typical of star-bursts. The 
resonance line 0[II]A3727 may also be detectable in stacked SDSS 



spectra since the host-galaxies of 



> 2A absorbers will 



be very close to the line-of-sight. The detectability of this emission 
line should be a strong function of VK/*^^®® but will also obviously 
depend on the host-galaxy metallicity. Secondly, the Nal D dou- 
blet ought to be detectable in the same stacked spectra. Thirdly, the 
colours of the stacked light distribution should be a strong function 
of Wr^'^^^". The preliminary results of S. Zibetti (private commu- 
nication) confirm this last prediction. 
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